Non-equilibrium fluctuations in a driven stochastic Lorentz gas 
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We study the stationary state of a one-dimensional kinetic model where a probe particle is driven 
by an external field £ and collides, elastically or inelastically, with a bath of particles at temperature 
T. We focus on the stationary distribution of the velocity of the particle, and of two estimates of 
the total entropy production As to t- One is the entropy production of the medium As m , which is 
equal to the energy exchanged with the scatterers, divided by a parameter 8, coinciding with the 
particle temperature at £ = 0. The other is the work W done by the external field, again rescaled 
by 9. At small £, a good collapse of the two distributions is found: in this case the two quantities 
also verify the Fluctuation Relation (FR), indicating that both are good approximations of Astot- 
Differently, for large values of £, the fluctuations of W violate the FR, while As m still verifies it. 

PACS numbers: 05.40.-a,02.50.Ey,05.70.Ln 

I. INTRODUCTION 

The program of statistical mechanics, consisting in deriving macroscopic properties of a system from the elementary 
interactions of its constituents, is far from being fulfilled in out-of-equilibrium conditions. In particular, the dissipation 
of energy in such cases prevents the use of general equilibrium results, and forces one to rely on a case-by-case model- 
dependent description. 

One of the few general results for systems far from equilibrium has been established in a series of works, starting from 
the seminal papers by Evans and coworkers [l|. It consists in a family of relations, generally referred to as Fluctuation 
Relations (FR), which are very similar in form, but concern different quantities (e.g. phase space contraction rate, 
entropy production, heat, work, etc.) and/or different dynamical regimes (e.g. transients, stationary states, etc.), as 
well as different kinds of non-equilibrium systems, either deterministic or stochastic (see [2| and references therein). 

In the framework of stochastic dynamics, the microscopic definition of the entropy production relies on the knowledge 
of the path probabilities for the model JjMj. The study of the fluctuations of this quantity plays a central role in the 
characterization of small systems [ToL Hlj . The connection with macroscopic quantities that are reasonably related 
to the thermodynamic concept of "entropy produced by the system" has to be carefully investigated in each specific 
model, as for instance in |12h14| . In particular, the simple recipe, useful near equilibrium, where the macroscopic 
entropy production is expressed as the work done by external forces divided by the temperature [15ft. is hardly of use 
far from equilibrium, often because it is not clear which parameter plays the role of temperature [16j. For instance, 
in granular gases, due to the dissipativc character of interactions, the kinetic temperature T g of the microscopic 
constituents does not always have a thermodynamic role [l7| . 

In order to address such issues, we investigate the dynamics of a stochastic Lorentz-like model where a particle 
is interacting with random scatterers and is subjected to an external force. Stochastic Lorentz models have been 
previously studied, focusing on transport properties, e.g. on normal or anomalous diffusion, for instance in (]~8l — 
[231 ] . Our model is based on the following ingredients: 1) the presence of an external field £ accelerating the probe 
particle, 2) scatterers of finite mass which are randomly and uniformly distributed in space and move with random 
velocities as extracted from a thermal bath at temperature T (it is therefore more reasonable to call them "bath 
particles"), 3) collisions which can also be inelastic (the system always reaches a stationary state), 4) a uniform 
collision probability which is inspired from the so-called Maxwell- molecules models [24j , and which helps in simplifying 
analytical calculations. Such ingredients provide a system in a non-equilibrium stationary state (NESS), due to the 
presence of a finite stationary current of particles. We do not consider any time-dependent experimental protocol 
neither any transformation between different NESS, that is, we study the system at fixed value of the external 
parameter (the external field). For each value of the external field, we start our measurements when the system has 
already reached the corresponding NESS. Therefore, in our case, there is not excess heat due to transient dynamics 
between different NESS and the total heat exchanged equals the power dissipated to sustain the NESS. 

We show that, in this model, the microscopic entropy production is well approximated, at large times, by the 
energy exchanged with the bath, divided by a "temperature" 6 which is that measured for zero field, £ = 0. Such 
a temperature turns out to be different from the temperature of scatterers T, unless the collisions are elastic, and 
coincides with the one of the probe particle in the unperturbed process, in agreement with what recently found in [25j : 
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at small values of the field, a macroscopically accessible quantity well approximates the entropy production, and that 
is the work done by the external field, divided by the same temperature 9. 



II. THE MODEL 



We consider an ensemble of probe particles of mass m endowed with scalar velocity v. Each probe particle only 
interacts with particles of mass M and velocity V extracted from an equilibrium bath at temperature T: such scatterers 
are distributed randomly and uniformly in space and can hit the particle only once. The last condition is important 
to guarantee unbounded motion and molecular chaos even in one dimension, and can be thought as the effect of 
bath particles moving in two (or more) dimensions, while the probe particle can only move along a one-dimensional 
track. Inspired from Maxwell-molecules models [24| and to their inelastic generalization [26U28I ] . we assume that the 
scattering probability does not depend on the relative velocity of colliders. Velocity of the particle changes from v to 
v' at each collision, according to the rule: 

v' = 7 v + (1 - i)V, (1) 

where 

7=^ (2) 

with £ = m/M, and a is the coefficient of restitution determining if the collision is elastic a = 1 or inelastic a S [0, 1). 
The velocity V of the bath particles is a random variable generated from a Gaussian distribution with zero mean and 
variance T/M: 



W) = ^exp(-^). (3) 

In addition, the probe particle is accelerated by a uniform force field m£. The resulting system can be assimilated to 
a Lorentz gas model, where free flights in external field are interrupted by random collisions with scatterers. 

The model is summarized by the linear Boltzmann equation for the evolution of the velocity distribution of the 
probe particle 

T c d t P(v,t)+T c £d v P(v,t) = 

-P(v, t) + j^— Jdu P(u, t)P s {^Y^) , (4) 

where t c is the mean collision time. In the following, we shall compare analytical predictions with numerical simulations 
of Eq. (U]). This equation, restricted to the particular case a = 1 and m = M (that is 7 = 0), has been recently 
studied in 1291. 



III. ENTROPY PRODUCTION 



One of the most interesting features peculiar to out-of-cquilibrium stationary dynamics is that a finite rate of 
entropy production can be measured in the system. Microscopically, such a quantity is related to the violation 
of detailed balance and gives a measure of how the probability of observing a forward trajectory differs from the 
probability of observing the time-reversed one [H, H| • From the macroscopic point of view, the entropy production is 
related to the presence of currents going through the system, due to spatial gradients or due to the action of external 
driving forces [HI, US]. In simple examples 0, Hl|, a bridge between the two points of view can be verified, where 
the microscopic entropy production turns out to be proportional to the product of a flux by a force. The constant 
prcfactor, in driven systems in contact with a reservoir, is often found to be the bath temperature. However this 
cannot be the general situation: for instance, for arbitrarily strong external fields, the temperature of the system (e.g. 
the kinetic one) may be far from that of the bath, or in extreme cases, cannot even be defined. 

The stochastic process considered here consists of two parts: a deterministic evolution, due to the action of the 
external field, plus a random contribution, due to the collisions with the scatterers. Here, the deterministic process 
does not contribute to the entropy production in the system, since the probability of a free fall for the particles is 
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symmetric under time-reversal. It is then convenient to rewrite Eq. (|4]) as a Master Equation where the transition 
rates describing the stochastic collisions between particles explicitly appear 



T c d t P(v,t)+T c £d v P(v,t) 

dv w(v\v )P(v , i) 



dv' 'w[v'\v)P{v, t), 



with 



w(v'\v) 



7 y/2KT/M 



exp 





v' + jv 


2 j 


I ~2T 


L 1 - 7 . 





(5) 



(6) 



To obtain this result we have put in Eq. Q the form (|3|) for the distribution of the scattering particles (see Ref [32[ 
for the general case, with different interaction kernels and arbitrary dimension). 

For the process described by Eq. ([5]), we can explicitly write the probability density P[{z)(i)}|t>o] of observing the 
trajectory {v(s)}l =to in the interval [to,t], with initial and final values v(to) = Vq and v(t) = v tl respectively. The 
total entropy production associated with each trajectory is defined as 



As tot (i) = log 



P[{v(t)}\v Q }P(v ) 



= As TO (i)+log 



P(-vtY 



(7) 



s )}s=* the time-reversed path, P(v) the stationary distribution and As m the entropy production 



with {v(s) = —vCt 

of the medium [33j, where as medium we refer to the ensemble of scatterers. Notice that in the reversed protocol 
the electric field does not change the sign, £ = £ . Along a trajectory where N c collisions occur at times tj (with 
j £ [1, Ac]) and velocities change from Vj = v(tj) to v'j = v(tj) one has 



nf^ 1 Pnc(A^>K> J ) 



t„ - t 



(8) 



Since in the Maxwell 



P[Mt)}\vo] Uf=lPnc(At 3 )w( 

where P nc (Atj) is the probability that no collision occurs in the time interval Atj 

model the collision probability is independent of the particle velocity, the time intervals Atj between successive 
collisions are distributed according to a Poissonian process with P nc (Atj) ~ exp(— Atj/r c ). In the ratio between 
probabilities appearing in Eq. ([8]). the contributions due to free flights in the backward trajectory exactly cancel those 
coming from the forward trajectory. It is interesting to mention that such cancellation is not a specific feature of the 
chosen uniform collision probability: it can be verified also for other kinds of interactions, e.g. for hard spheres. 

The stochastic entropy production of the medium is then expressed in terms of single event contributions 6s m {v', v), 
where the particle changes its velocity from v, before the collision, to v', after the collision, 



JV e (t) 

As m (t) = 5s m {v'j,Vj). 
Here N c (t) is the number of collisions occurred up to time t and, from Eq. ([5]), 



(9) 



SSmiv'pVj) = log 



W(-Vj \ - it) 



M(l + 7) , ,2 2, 

SEcouty'pVj) 



where 8E co u(v'p Vj) = m/2(v' 1 2 — in) is the energy received from the scatterers in a collision and 



9 = T( 



1-7 



1 + a 



1 + 7 ' 1 + 2C - a 

is the kinetic temperature of the probe particle in zero field, as demonstrated below (see Eq. 



(10) 



(11) 
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For the average rate of entropy production in a stationary state we can write 

(s m (t)) = lim (As 7 (t)) = - (6s m (v',v)) 

1 M(l+ 7 ) , 2 2 

-({V )post ~ {V )), 



r c 2T(1 - 7 ) 

(12) 

where (• • • ) pos t denotes the average over the distribution of post-collision velocities. In order to compute explicitly 
the average quantities appearing in Eq. (|12p . we need to solve the Boltzmann equation (U) in the stationary limit. 



IV. STATIONARY SOLUTION 

Defining the Fourier transform of the probability density P(k,t) = dve lkv P(v,t) and similarly for Pg(fc), the 
integral equation Q in Fourier space in the stationary limit assumes the convenient form: 

- ik£r c P(k) = -P{k) + P( 7 fc)Ps((l - 7 )fc). (13) 

Let us start by considering the simplest case 7 = 0. This corresponds to ( = a and it has been considered in [29l |34|] . 
for elastic particles with m = M . Since P(0) — 1, from Eq. ([lU|) we have 

Hk) = ^r-TPs(k) = -4^T exp (-^jkA . (14) 



1 - i£r c k w 1 - i£r c k 1 \ 2M 
The inversion of the Fourier transform gives the sought distribution under the form of a convolution: 



p(^ v ) — bJ— I c?ucxp(— a(v — u) 2 ) x exp(— bu), (15) 



with b = and a = ^p. The result of the convolution is that the tails for large v are exponential and for negative v 
are Gaussian. In the limit of infinite mass of the scatterers, M — >• oo, the distribution becomes one sided and reads: 

P{v)=6{v)^-exp(-^]. (16) 



t c £ \ t c £ 

In Fig.[T]thc analytical prediction of Eq. (|15p is compared with the pdf obtained in numerical simulations (see below). 
For arbitrary values of 7, equation (| 13[) can be easily solved in the case of zero field £ = 0, where one has 

P (k) = exp ( -^-k 2 ) , (17) 
\ 2m / 

with a temperature 9. Notice that Pq(v) is a Maxwcllian, but with a "temperature" 6 which differs from T; one has 
6 = T only when a = 1. In the presence of non-zero field, no analogue of the closed formula (|15j) can be written. 
However, in this case one has access to all the moments of the distribution. Indeed, assuming analyticity around 
k = 0, for small enough k one can write the expansion 

p(fe) = f;^M„ ) as) 

^-^ TV. 

where fi n = (v n ). Upon substituting expression (TT5|) in Eq. (fT3")) . and equating equal powers of fc, all the moments of 
the distribution can be obtained. In particular, using Eq. ([3]), we have 

Mi = -p—— (19) 
1-7 



(20) 
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FIG. 1: (Color online) Comparison between the theoretical prediction of Eq. (|15[1 (solid lines) and the numerical simulations 
(symbols) for the stationary distribution in the case 7 = 0, r c = 1, C = 0.5, T = 1 and for different values of the external field 
£ = 0.5,2,8. 



It is interesting to note that the conductibility (v)/£ = r c /(l — 7) = r c (l + C)/(l + a ) increases when the system 
becomes more inelastic (a reduced) or when the mass M of the bath particles is reduced. 

We are now ready to compute the average entropy production. Indeed, substituting Eq. ([20)1 in Eq. (|12p , and using 
the fact that the distribution of post-collisional velocities is by definition P pos t(k) = P(7fe)JPg((l — 7)&), we can write 



.(*)> 



(1+7) M 
(1-7) 2 T 



t c £ 2 > 0. 



(21) 



Expression (|2ip can be related to the macroscopic quantities present in the system, namely the external field £ and 
the current velocity (v). Indeed, using Eq. (JTTJ, Eq. (|2"Tj) can be rewritten as 



m£(v) 



> 0. 



(22) 



(1-7) 

Now let us consider the average work done by the external field along a trajectory that spans the time interval [0,t], 
W(t) = f* Fv(s)ds, with F = m£: 

F(v) 



y 1 (W(t)) 
lim — 

t->oo t 



»(*)), 



(23) 



i.e. the average macroscopic work of the field divided by the temperature 8 corresponds to the average entropy 
production. Notice that in Eq. (|23|) the "right" temperature is neither the bath temperature T, nor the kinetic 
temperature of the probe particle in the presence of the field 



T g = m(fi 2 - Mi 



(Tc£? 

'(1-7 2 V 



(24) 



but it is that of the unperturbed, not accelerated system. The quantity 9 represents an energy scale in the system, 
which depends on the several parameters defining the model, namely temperature of the scatterers T, mass ratio £ 
and restitution coefficient a. It is equal to the kinetic temperature of the particle only in the absence of external field. 
In this case, if the interactions are elastic, it equals the bath temperature T, as expected. In general, for non zero 
field, the relation between 9 and the mean square velocity of the particle is expressed by Eq. (l24l) . 



NUMERICAL SIMULATIONS 



In order to obtain numerically the stationary distribution P(v) which solves Eq. (U), we simulate the dynamics of a 
single particle subject to a constant acceleration £ and to inelastic collisions with the scatterers, and average over 10 4 
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FIG. 2: (Color online) PDF of entropy production, work and rescaled work. Parameters: 7 = 0, t c = 1, £ = 0.5, T = 1, £ = 0.5 
and t = 16t c . The continuous line shows the Gaussian fit of p(As m ). 



realizations, for general 7. Time is discretized in intervals 5t, for a total duration 8tN t . The particle is accelerated 
for St under the influence of the constant field £ and then a collision with a scatterer is realized with probability 
Pcoii — St/r c . The distinguishing feature of Maxwell molecules [24| is that the probability of colliding is independent 
of the velocity of the particle itself. Then, the particle is again uniformly accelerated by the field £ and the whole 
procedure is repeated. This ensures an average collision time r c , which is set as a parameter of the simulation. For 
each collision the velocity of the scatterers are extracted from a constant Gaussian distribution, Eq. <j3j> , and the 
collision with the probe particle is realized with the inelastic rule, Eq. (JXJ) . 

In the numerical simulations we also studied the diffusive properties of the system to find that the mean square 
displacement with respect to the average motion, e?2 = ([x(t) — (x(t))} 2 ), displays a ballistic behavior d,2 ~ t 2 on short 
time scales with a crossover to a diffusional dynamics d? ~ t at large times. Such a ballistic-to-diffusive scenario, 
not intuitive in the presence of an accelerating field (l9l.]22|. is consistent with the autocorrelation of the particle's 
velocity which, up to numerical precision, decays as a single exponential. 



VI. THE ROLE OF INELASTICITY 



A peculiarity of this model is that the inelastic collision rule can be always mapped onto an clastic one: indeed we 
can always set a = 1 and change the mass ratio £ = m/M in order to keep constant the parameter 7 which enters the 
collision rule. In practice, if also m is kept fixed, we are changing only M and thus the width of the distribution in 
Eq. ([3]). Doing so we find no relevant qualitative change in the physics of the system. In particular, let us notice that 
the entropy production in Eq. (|21[) depends only on 7 and vanishes with the field £ , also with inelastic collisions. 

The fact that, for the present model, the entropy production is not affected by the inelasticity of collisions is 



in agreement with the findings of 35|, where the dynamics of a single massive intruder in a diluted granular gas 



of inelastic particles satisfying the molecular chaos hypothesis was studied. In that case the dynamics of the probe 
particle was well described by a single linear Langevin equation, so that the entropy production was zero by definition. 
Without passing through a Langevin description, also the Maxwell model presented here assumes the molecular chaos 
hypothesis for the surrounding sea of scatterers. In both cases each collision of the probe particle is performed with 
a velocity extracted from a constant distribution independent of the previous history of the system: despite the 
inelasticity of collisions, reversibility is guaranteed by the molecular chaos hypothesis for the medium. As it is clear 
from Eq. (|21|). in our model all irreversible effects are generated by the external field. 

There is only one special case in which the entropy production allows us to distinguish the inelastic from the elastic 
interaction. This is the case of infinitely massive scatterers, namely M — > 00, which implies C = and 7 = —a. 
Keeping finite the width of the velocity distribution of scatterers, T/M = \ tne average entropy production then 
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FIG. 3: (Color online) The symmetry relation (|26[) obeyed by the entropy production As m (t) and by the rescaled work W(t)/0 
obtained in numerical simulations with parameters: 7 = 0, t c = 1, C = 0.5, T = 1, £ = 0.5 and t — 16r c . Inset: collapse of the 
rate functions of the distribution of As m (t) for large times t = 128r c (black circles), t — 512r c (red squares) and t = 1024r c 
(green diamond). 



reads 



.(*)> 



1 - a t c £ 2 



(25) 



namely it vanishes in the elastic case. This case is quite peculiar because its stationary P(y) has finite (v) but infinite 
(v 2 ) [291 ] . For any other parameters the model has finite or zero current and finite energy. 



VII. FLUCTUATIONS OF ENTROPY PRODUCTION AND WORK 



By definition the total entropy production ([7| must satisfy the FR for any value of t: 

p(As to t(t) = x) 



log 



p(As to t(t) = -x) 



(26) 



with p(Astot(t) = x) the probability density that in a time interval [0,t] the entropy produced by the system is 
x. At large times t, one usually neglects the term \og[P(va) / P(—vt)] which - in the definition of Astot(t) - gives a 
contribution of order 0(1), and looks for the fulfillment of Eq. (f2l)]) using As tot (i) w As m (t), which should be the 
leading part of order 0{t). Our numerical simulations suggest that this is the case also in this model, for any value 
of the field £ . As seen in Eqs. (|9|)- (fT0|) . As m (t) coincides with the energy that the probe particle loses when colliding 
with the bath particles, divided by the kinetic temperature of the particle itself measured at zero field. In Fig. [2] the 
distribution of As m (t) is plotted for a certain value of t: it can be clearly seen from the superimposed fit a remarkable 
deviation from gaussianity. 

Inspired from the equality of the average value of the entropy and of the work produced along a trajectory, see 
Eq. (|23|) . wc measured in numerical simulation also the fluctuations of the work done by the external field on the 
particle, W(t). A remarkable finding is that, for small values of £, the PDF of the work can be collapsed on the PDF 
of the entropy production by exploiting the "unperturbed" temperature as a scaling parameter. In such cases, the 
measure of work fluctuations provides an alternative way to measure the temperature of the unperturbed system: 
indeed, as also verified in Fig. El when W(t)/6 satisfies the FR, one also has \og[p(W(t) = x)/p(W(t) = —x)] = x/9. 

By increasing the value of £ while keeping fixed all the other parameters of the simulation we find the following 
differences with the small £ situation: a) the stationary PDF of velocity has a larger exponential tail and a very 
asymmetric shape; b) the PDF of As m (t) also becomes very asymmetric and close to has a shape very different 
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FIG. 4: (Color online) PDF of the quantities As m , W/6 and AE/6 measured in numerical simulations with parameters: 7 = 0, 
r c = 1, £ = 0.5, T = 1, £ = 8 and t = 16r c . Inset: zoom of the region around x = for the probability distribution of As m . 



from a Gaussian one, see inset of Fig.|4] c) the PDF of W(t) cannot be collapsed on the PDF of the entropy production 
with a simple rescaling; d) As m (t) still verifies the FR, while W(t)/0 does not. 

The breaking of the FR symmetry by the fluctuations of work at large values of £ can be interpreted considering 
the balance equations of the energy absorbed and lost by the probe particle within a time window. In particular, 
we have that W{t) + AE coU (t) = AE(t), with AE(t) = f[v 2 (t) - v 2 (0)] and AE coU (t) = J^f^ ^colliv'^Vj) and 
therefore, exploiting the relation (jlOp . we can write: 

As m (i) = \ [W(t) - AE(t)] , (27) 

which puts in evidence that discrepancies between the distribution of As m and W/6 are due to the "boundary term" 
AE. Although AE(t) is of order C(l), it is known to be dangerous for FR, even at large times, when its distribution 
has exponential (or larger) tails [36|, 133], as in our case (see Fig. 2]). 

In summary, we have discussed a non-equilibrium kinetic model, simple enough to let most of the calculations 
accessible, such as the entropy production or the moments of the stationary distribution, but still displaying interesting 
properties, for instance strong non-Gaussian behavior and non-trivial dependence on the external field. In particular 
we have seen that the Fluctuation Relation is valid for the distribution of the energy AE co u lost in collisions with 
the bath, and it also holds for the work done by the field, W(t), provided that £ is low enough. For the fulfillment of 
the FR, both quantities have to be divided by a temperature 9 which is different from the bath and the particle ones, 
but coincides with the latter at £ = 0. We remark that such an observation is highly non-trivial: it suggests that 
extreme care must be used when entropy production is defined on phenomenological grounds, where one is usually 
tempted to use more "reasonable" temperatures (e.g. the bath or the system ones), forgetting the complexity of 
far-from-equilibrium systems. 
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